Single Stranded Breaks Relax Intrinsic Curvature in DNA 
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The macroscopic curvature of double helical DNA induced by regularly repeated adenine tracts 
(A-tracts) is well-known but still puzzling. Its physical origin remains controversial even though 
it is perhaps the best-documented sequence modulation of DNA structure. The recently proposed 
compressed backbone theory (CBT) suggested that the bending can result from a geometric mis- 
match between the specific backbone length and optimal base stacking orientations in B-DNA. It 
predicted that the curvature in A-tract repeats can be relaxed by introducing single stranded breaks 
(nicks). This effect have not been tested earlier and it would not be accounted for by alternative 
models of DNA bending. This paper checks the above prediction in a combined theoretical and 
experimental investigation. A series of nicked DNA fragments was constructed from two mother 
sequences with identical base pair composition, one including A-tract repeats and the other being 
random. The curvature was tested experimentally by gel mobility assays and, simultaneously, by 
free molecular dynamics (MD) simulations. Single stranded breaks produce virtually no effect upon 
the gel mobility of the random sequence DNA. In contrast, for nicked A-tract fragments, a regular 
modulation of curvature is observed depending upon the position of the strand break with respect to 
the overall bend. As predicted by CBT, the curvature is reduced, with maximal relaxation observed 
when the nick occurs inside A-tracts. These experimental results are partially reproduced in MD 
simulations of nicked A-tract repeats. Analysis of computed curved DNA conformations reveals 
modulations of local backbone length as measured by distances between some sugar atoms, with 
maximal backbone compression in the narrowing of the minor groove inside A-tracts. The results 
lend additional support to CBT versus alternative mechanisms of intrinsic curvature in DNA. 



I. INTRODUCTION 



The intrinsic curvature in DNA was discovered about 
twenty years ago for regular repeats of A n T m , with n + 
m > 3, called A-tract o 1 ' 2 ' 3 . Several profound reviews ex- 
ist of the large volume of related data accumulated dur- 
ing the last decades^aii£iZi2i&. Due to its outstanding role 
in genome functioning and a very controversial experi- 
mental behavior, the sequence dependent DNA bending 
continues to attract great attention. The hottest point of 
debates is the molecular mechanism of this phenomenon. 
Already in the seventies, it was realized that the DNA 
double helix is not a rigid cylinder, but is bendable de- 
pending upon direction and the base pair sequenc o 10 ^ 1 . 
By postulating that stacking in ApA steps is intrinsi- 
cally non-parallel the wedge model predicted that A- 
tracts repeated in phase with the helical screw should 
bend DNAJ^. This proved true, but the model later ap- 
peared unsatisfactory. It turned out that X-ray DNA 
structures are often bent in random sequences whereas A- 
tracts are exceptionally straigh t ^Wi? . With the grow- 
ing number of sequences checked, non-zero wedges had 
to be introduced for all base pair steps 1 ^ 7 *^, never- 
theless, experimental counter-examples were eventually 
found where bending could not result from accumula- 
tion of wedgesA&SH. Several theories tried to reveal other 
physical factors involved in DNA bending along with base 
stacking. The most popular junction model proposed by 
the Crother's group2*2i originated from an idea that a 



bend should occur when two different DNA forms are 
stacked 2 ^. If poly (d A) .poly (dT) double helix had a spe- 
cial B' form as suggested by some data^ 3 - the helical axis 
should be kinked when an A-tract is interrupted by a ran- 
dom sequence. Drew and TraversS 4 *^ apparently were 
the first to notice that narrowing of both DNA grooves 
at the inner edge of a bend is a necessary and suffi- 
cient condition of bending. They, and later Burkhoff and 
Tullius 2 ^, considered the preference of narrow and wide 
minor groove profiles by certain sequences as the possible 
original cause of this effect. Similar ideas were considered 
within the context of the junction modelSl. In addition, 
the possible involvement of solvent counterions in DNA 
bending was discussed long agoSS, and it has been exten- 
sively studied in the recent years^ffiiSi. Detailed analy- 
ses of all these theories including comparison with exper- 
imental data can be found in the recent literatur oW' 32 . 
Here we note only that, in spite of all efforts, the contro- 
versies remain and continue to accumulat o 33 i 34 . 

Conformational calculations were always intensively 
used for testing various hypothesis concerning DNA 
bendin g 11 ! 21 ^ 5 ^^ 37 ! 3 ^ 3 ^ 40 . The progress in method- 
ology of molecular dynamics (MD) simulations of nu- 
cleic acids^i recently made possible realistic modeling 
of A-tract DNA fragmenta 27 i 42 i 43 i 44 . In straightfor- 
ward simulations of bending, the initial straight sym- 
metrical DNA is allowed to bend spontaneously with 
no extra forces applied, that is due to generic atom- 
atom interactions 2 ^ 3 ^ 4 ^ 4 ^ 4 ^. Such computational ex- 
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periments are very demanding because they involve rela- 
tively long DNA fragments and require multi-nanosecond 
(at least) trajectories to reveal reliable trends toward 
bent states. Nevertheless, for preselected A-tract se- 
quences, it was possible to demonstrate spontaneous 
emergence of intrinsic curvature in the course of sim- 
ulations as well as attracting character of bent states, 
which is necessary for static curvature. These simula- 
tions resulted in a new hypothesis of the physical origin 
of intrinsic curvature which we refer to as the compressed 
backbone theory (CBT) 3 WF. 

In a certain sense, the CBT continues the line of ear- 
lier views that attributed the origin of bending to the 
properties of the DNA backbon o 24 i 25 i 26 i 27 . It postulates 
that the intrinsic curvature in B-DNA results from a 
geometric mismatch between the length of the sugar- 
phosphate backbone and the base pair stacking. Starting 
from the first published drawing of the double helis^, 
its two backbone strands are viewed as regular spiral 
traces that wind along the surface of the cylindrical B- 
DNA core formed by stacked base pairs. Because reg- 
ular spirals represent the shortest lines that join two 
points on a cylindrical surface, and because the back- 
bone strands are covalently linked to bases, an ideal B- 
DNA can exist only if the average backbone length ei- 
ther matches exactly or is shorter than that dictated by 
the base pair stacking 4 ^. The CBT assumes the oppo- 
site. It postulates that, under physiological conditions, 
the equilibrium specific length of the B-DNA backbone, 
considered as a restrained polymer attached to a cylin- 
drical surface, is slightly longer that in the canonical B- 
form. In physics, similar mismatching relationship is of- 
ten called geometric frustration 4 ^. The backbone tries 
to expand and "pushes" stacked bases while the stack- 
ing interactions oppose this. Eventually, the backbone 
finds a "compromise" by deviating from its regular spiral 
trace, which causes quasi-sinusoidal modulations of the 
DNA grooves. Concomitant base stacking perturbations 
cause local bends that accumulate to macroscopic static 
curvature when the period of these modulations corre- 
sponds to an integral number of helical turns. The CBT 
considers B-DNA double helices with any sequences and 
it may have a number of interesting consequences. It 
suggests, for example, that, even under normal temper- 
ature, the B-DNA double helix can occur in a glass-like 
state characterized by very long relaxation timed. The 
latter prediction agrees with some recent experimental 
d&taHtSML. Under certain assumptions, one can even 
imagine the possibility of a glass-like aging of a single 
DNA molecule in vivo, with direct relation to the prob- 
lems of biological aging. 

One of the possible experimental tests of CBT is 
based upon the predicted effect of single-stranded breaks 
(nicks) upon the intrinsic curvature^. Reversible single 
chain breaks in DNA occur during diverse biochemical 
processes, including general and site-specific recombina- 
tion, replication, and DNA repair. Therefore, structural 
perturbations due to nicks and single-stranded "gaps" 



have been earlier characterized in a number experimen- 
tal i -| lj r, ' .-.r, r„. .-,7 r,s .-,1. . ,u r,r, «7 ,<» Ac _ 

cording to all these data, nicks perturb the free DNA 
very little. The corresponding X-ray structured 
are straight and similar to those of intact duplexes 
with analogous sequences. This agrees with early 
physicochemical test a 52 ' 5 !?' 54 as well as NMR studies in 
solution 56 i 57 i 65 i 66 i 67 . The main detectable effect of a 
nick is local melting or fraying, especially at low ionic 
strength, mild denaturating conditions, or elevated tem- 
peratures. This increases the isotropic flexibility of DNA 
and reduces its gel mobility£2i£ii£2i££. In normal condi- 
tions the nick fraying is small and it almost undetectable 
below 10 o CSMa. 

Based upon these experimental data and the earlier 
views of the origin of intrinsic curvature one should ex- 
pect that nicks in phased A-tract sequences can either 
increase bending or have zero effect. The increase may 
be expected if the curvature is caused by external forces, 
as in electrostatic counterion models, because the flex- 
ible hinge at a nick can allow these forces to increase 
the bend angle in their direction. A zero effect is ex- 
pected for the wedge and junction models because nicks 
arguably do not affect the A-tract structured and do 
not cause wedges detectable in gel migration assayd. 
In contrast, if the intrinsic curvature is really caused by 
the backbone compression, single-stranded breaks should 
relax it, which should increase the gel mobility of curved 
DNA fragments. This effect is opposite to that of the 
increased isotropic flexibility, moreover, the latter can be 
eliminated by reducing the temperature so that these two 
phenomena cannot be confused. In addition, because the 
degree of the compression should vary regularly between 
the inner and outer edges of the bent DNA, relaxation 
produced by a single nick is expected to change systemat- 
ically according to its position with respect to the overall 
bend. 

Here we present the results of experimental and simu- 
lation studies carried out according to the above plan. 
Two 35-mer sequences from our recent comparative 
investigation^ 2 , were used for constructing two series of 
nicked DNA fragments. The first sequence contained an 
A-tract repeat motif that was preselected in MD simula- 
tions to reproducibly induce strong static curvature. The 
second sequence was "random", but with the same base 
pair composition. Backbone breaks were introduced in 
one strand at different portions to span approximately 
one helical turn. The curvature was probed by PAGE as 
well as long-time MD simulations. In experiment, a no- 
ticeable increase in the gel mobility of the nicked A-tract 
repeat is observed depending upon the nick position. In 
contrast, no such effect is found for the random frag- 
ment. MD simulations fail to reproduce experimental 
effects qualitatively, but qualitatively they also confirm 
that single stranded breaks interfere with the A-tract in- 
duced curvature. The results qualitatively agree with 
CBT, but they cannot be accounted for by other models. 
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II. MATERIALS AND METHODS 

A. Oligonucleotides and construction of 5'-labeled 
DNA probes 

The double stranded DNA containing four A-tracts 
flanked poly (d A) .poly (dT) termini were constructed by 
annealing the synthetic oligonucleotide B (the top line in 
Table Pi with its complementary (bottom) strand. Dou- 
ble stranded DNA fragments containing single stranded 
breaks (nicks) were constructed by annealing two shorter 
oligonucleotides (Table HJ with the same bottom strand. 
In the same way a series of nicked duplexes was con- 
structed from the reference random sequence fragment 
S, with its base pair composition identical to that of B. 
The sequences of the oligonucleotides used for the con- 
struction of the DNA fragments are assembled in Table 
U The fragment codes in this table use the following 
mnemonics. The capital 'B' and 'S' stand for "bent" and 
"straight". The small letter prefix 'n' stands for "nick". 
The numbers indicate the shift of the nick position with 
respect to the center of the sequence. 

In all the constructs, one of the oligonucleotides an- 
nealed was labeled with T4 polynucleotide kinase and [ 32 - 
P]-ATP. In most cases the label was attached to the 5'~ 
end of the bottom (continuous) oligonucleotide whereas 
the two strands at the nick position were terminated 
with 3'OH and 5'OH groups. For some sequences, 5'- 
phosphorylated nicks were additionally constructed to 
check the effect of phosphate charges. In this case the 
continuous strand carried 5'-OH while [ 32 -P] was intro- 
duced in the corresponding shorter partner. The anneal- 
ing was carried out by incubating 300 nM of the unla- 
beled oligonucleotide (s) with 30 nM of end-labeled ones 
for 3 min at 80° C in 20 mM Tris-HCl (pH 8.0), 150 mM 
NaCl, and then allowing them to cool slowly. To dis- 
tinguish between the two nick series, in the text, we add 
small letter suffixes 'o' and 'p' to the codes given in Table 
Q] Thus the whole two nick series of the bent DNA are 
referred to as nBo and nBp, respectively, while the in- 
dividual nicked fragments are denoted as nB-2o, nB-2p, 
etc. 



B. Gel mobility assays 

The mobility of the DNA fragments was analyzed in 
12% gels (acrylamide to bis-acrylamide, 29:1) buffered 
with 90 mM Tris-borate, 1 mM EDTA, pH 8.6. In order 
to check the effect of Mg 2+ ions the buffer was supple- 
mented with 10 mM MgCl 2 in the absence of EDTA. Gels 
were pre-run under constant power until stabilization of 
the current. Labeled DNA in a buffer containing 20 mM 
Tris-HCl (pH 8.0), 50 mM NaCl, 4% Phicoll-400 and xy- 
lencianol was loaded onto the gel. The electrophoresis 
was performed under constant voltage and constant tem- 
perature of 4°C. The dried gels were exposed to storage 
phosphor screens and visualized on a 400S Phosphorlm- 



ager (Molecular Dynamics). 



C. Calculations 

MD simulations were carried out for the central 35 base 
pairs of six nicked DNA fragments in Table Q] constructed 
from the A-tract repeat sequence (nB-4o, nB+6o). 
The corresponding calculations for the two mother frag- 
ments without nicks (B and S) were reported earlier* 3 ^. 
For consistency with the previous studie o^ 2 i 4, ^ 4 ff , we used 
the AMBER98£2iia force field and TIP3P watenii. Tra- 
jectories generally started from the fiber canonical B- 
DNA conformation 7 ^ and were continued to 20 ns. In 
a few cases, dynamics was also run starting from the 
canonical A-DNA conformation. In total, these simula- 
tions sampled from more than 0.2 ms of all-atom dynam- 
ics for nicked and intact 35-mer DNA fragments with the 
same A-tract repeat sequence. 

Molecular dynamics simulations were carried out with 
the internal coordinate molecular dynamics (ICMD) 
methodZ&Si adapted for DN A 75 i 76 with the time step of 
0.01 ps. In this approach, the DNA molecule has all 
bond length and almost all bond angles fixed at their 
standard values. The only variable bond angles are those 
centered at the sugar C1 ! ,...,C4 ! , and 04' atoms, which 
assures the flexibility of the furanose rings. In contrast, 
bases, thymine methyls, and phosphate groups move as 
articulated rigid bodies, with only rotations around sin- 
gle bonds allowed. The highest frequencies in thus ob- 
tained models are additionally balanced by increasing ro- 
tational inertia of the lightest rigid bodies as described 
earlies 7 ^?.. The possible physical effects of the above 
modifications have been discussed elsewher o i 78 . 

The so-called "minimal model" of B-DNA was 
use d 7 fti7ff . It includes only a partial hydration shell and 
treats counterion and long range solvation effects implic- 
itly by reducing phosphate charges to -0.5 and applying 
linear scaling of Coulomb forces. This model produces B- 
DNA structures very close to experimental data and has 
no other bias towards bent or non-bent conformations 
except the base pair sequence. The interest to implicit 
solvation in DNA simulations is long standing and it con- 
tinues in the literature^. Advantages as well as limita- 
tions of such approaches have been recently reviewed^. 

The minimal model is not meant to be generally ap- 
plicable and we have chosen it for the present studies 
because of the following considerations. We assume that 
effects it ignores, like sequence specific counterion bind- 
ing, play some role, but are not critical for the A-tract 
induced curvature. We have shown earlier that, in these 
conditions, the A-tract repeat motif we are using has a 
bent state with very strong attracting properties that al- 
low one to observe spontaneous transitions to stable bent 
conformations in the course of MD 3 ^. This property 
is rather exceptional and only due to it one could hope 
to obtain useful information from a comparison between 
nicked and intact DNA in MD. We also carried out sev- 
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TABLE I: Construction of DNA fragments. A-tracts are boldfaced and nick positions are marked by double quotes. 



Code 


Sequence 


B 


Ai 8 -AAAATAGGCTATTTTAGGCTATTTTAGGCTATTTT-Ti8 


nB-9 


Ai8-AAAATAGG"CTATTTTAGGCTATTTTAGGCTATTTT-Ti 8 


nB-6 


Ai8-AAAATAGGCTA"TTTTAGGCTATTTTAGGCTATTTT-Tis 


nB-4 


Ai8-AAAATAGGCTATT"TTAGGCTATTTTAGGCTATTTT-Ti 8 


nB-2 


Ai8-AAAATAGGCTATTTT"AGGCTATTTTAGGCTATTTT-Ti 8 


nB+0 


Ai 8 -AAAATAGGCTATTTTAG"GCTATTTTAGGCTATTTT-Ti8 


nB+2 


Ais-AAAATAGGCTATTTTAGGCTATTTTAGGCTATTTT-Tis 


nB+4 


Ai8-AAAATAGGCTATTTTAGGCTA"TTTTAGGCTATTTT-Tis 


nB+6 


Ai8-AAAATAGGCTATTTTAGGCTATT"TTAGGCTATTTT-Ti 8 


nB+11 


Ai8-AAAATAGGCTATTTTAGGCTATTTTAGG"CTATTTT-Ti 8 


S 


Ais-TTAGATAGTATGACTATCTATGATCATGTATGATA-Tis 


nS-9 


Ai8-TTAGATAG"TATGACTATCTATGATCATGTATGATA-Ti 8 


nS-6 


Ai8-TTAGATAGTAT"GACTATCTATGATCATGTATGATA-Ti8 


nS-4 


Ai8-TTAGATAGTATGA"CTATCTATGATCATGTATGATA-Tis 


nS-2 


Ai8-TTAGATAGTATGACT"ATCTATGATCATGTATGATA-Ti 8 


nS+0 


Ai8-TTAGATAGTATGACTAT"CTATGATCATGTATGATA-Ti 8 


nS+2 


Ai8-TTAGATAGTATGACTATCT"ATGATCATGTATGATA-Ti 8 


nS+4 


Ai8-TTAGATAGTATGACTATCTAT"GATCATGTATGATA-Ti 8 


nS+6 


Ai8-TTAGATAGTATGACTATCTATGA"TCATGTATGATA-Ti 8 


nS+11 


Ai8-TTAGATAGTATGACTATCTATGATCATG"TATGATA-Ti 8 



eral simulations for the mother 35-mer and 25-mer A- 
tract fragments by the Particle-Mesh Ewald method^! 
with full solvation and periodical boundaries. The cor- 
responding trajectories were continued to 5-10 ns, but 
they showed much less demonstrative bending dynamics. 
The origin of this difficulty is currently unclear, but it 
agrees with recent reports by some other authors^. Fi- 
nally, the minimal model provides for substantial savings 
in computations. The saving factor is critical because 
DNA bending is detectable only in relatively long frag- 
ments and it is likely to involve very slow motions. In 
spite of a few reported analogous simulations using full- 
scale solvatio n 27 i 42 such calculations would be too costly 
for the volume of sampling reached here. 



ally base pairs of the A-tract repeat. We preferred this 
randomized sequence to commonly used GC-rich straight 
fragments in order to keep the base pair content identical 
and reduce the noise that could cause small variations in 
gel mobility. In order to amplify the PAGE resolution of 
similarly bent nicked DNA, the 35-mer fragments were 
extended to 71 bp by adding poly (dA). poly (dT) tails. 
The tails continue the mother 35-mer A-tract repeat se- 
quence smoothly, which should reduce the possible per- 
turbations that could affect the comparison of PAGE re- 
sults with MD simulations that involved only the central 
35-mer nicked fragments. 



B. Relative PAGE Mobilities 



III. RESULTS AND DISCUSSION 
A. Construction of DNA Fragments 

The A-tract motif AAAATAG originally attracted our 
attention in MD simulations of the natural DNA frag- 
ment taken from the first curved DNA locus studied in 
vitr(&2&. The central 35-mer A-tract repeat for the bent 
fragments in Table Q] was constructed by repeating this 
motif four times and it had to be inverted to make the 
two DNA termini symmetrical. Such inversion should 
not affect bending^ but is essential for simulations be- 
cause the 3'- and 5 '-end A-tracts may represent qual- 
itatively different boundaries. In repeated simulations 
with this and similar A-tract fragments, the static cur- 
vature emerged spontaneously and it became more ev- 
ident as the chain length increasec~2ii£ . To obtain a 
reference non-A-tract DNA, we have re-shuffled manu- 



A representative PAGE plate of the random sequence 
series of non-phosphorylated nicks is shown in Fig. 
All nicked DNA exhibit similar mobilities close to that 
of the reference straight fragment. A weak retardation 
effect distinguishable for some nick positions is at the 
limit of experimental accuracy. Qualitatively similar pat- 
terns, but with stronger retardation were earlier reported 
for nicked DNA fragments under elevated temperature 
and/or mild denaturating condition oW'^ 1 '^ . We suppose, 
therefore, that some reduction of mobility in Fig. Q re- 
sults from the isotropic flexibility at nick positions that 
is strongly reduced at 4°C. This small retardation should 
be taken into account in the interpretation of other tests 
below. 

Fig. [2] shows the results obtained in the same condi- 
tions for the A-tract nicks. The pattern exhibited here 
is evidently more complex than that in Fig. On the 
one hand, for the central nick positions, the mobilities 
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FIG. 2: PAGE analysis of comparative mobilities of nicked 
DNA constructed from the A-tract repeat fragment. The con- 
tents of the lanes is as follows (see Table Q] and Methods for 
fragment codes). The lanes marked 'B' and 'R' ("reference") 
contain the intact B and its mixture with S, respectively. The 
series of lanes '-6', '+6' contain nB-6o, nB+6o, respec- 
tively. 
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PIG. 1: PAGE analysis of comparative mobilities of nicked 
DNA constructed from the random fragment. The contents of 
the lanes is as follows (see Table and Methods for fragment 
codes). The lanes marked 'R' ("reference") contain a mixture 
of intact duplexes S and B. The series of lanes '-6', '+6' 
contain nS-6o, nS+6o, respectively. 



are close to that of the mother fragment, which is qual- 
itatively similar to Fig. However, when the nick is 
moved farther from the center, the PAGE mobility grows 
and becomes higher than that of the mother curved frag- 
ment. It reaches maximum values for nB-4o and nB+6o, 
but for nB-6o it is reduced again. As a result, the mobil- 
ities of nB-6o and nB-4o appear close to those of nB+4o 
and nB+4o, respectively, and the overall pattern of nick 
bands in Fig. appears skewed-sinusoidal, with the pe- 
riod corresponding to that of the double helix. The in- 
crease of PAGE mobility indicates reduced DNA curva- 
ture and this effect can well be due to the relaxation pre- 
dicted by CBT. The amplitude of the reduction is very 
significant since it exceeds 30% of the curvature differ- 
ence between the S and B fragments produced by four 
A-tracts. This mobility profile suggests that the back- 
bone compression is maximal within A-tracts and mini- 
mal between them. Although this specific phasing is right 
opposite to our first guess^ it is accounted for by CBT 
and not by other theories as we discuss further below. 
However, there are a few alternative interpretations that 



should be considered first. 

The first evident factor that should be checked is the 
perturbed balance of phosphate charges. The two nick 
series displayed in Fig. Q]and Fig. [2 have one phosphate 
group less than the reference fragments B and S. Con- 
sequently, their mobilities may be uniformly reduced by 
around 0.7%, which would mean that the apparent reduc- 
tion of curvature in the nicked A-tract fragments nB-4o 
and nB+6o is even stronger than it is seems from Fig. |2 
This effect, however, is not distinguishable in Fig. Q]sug- 
gesting that it is too small and may be safely neglected. 
The influence of phosphate charges can be more subtle, 
however. The electrostatic models of DNA bending as- 
sume that it is the broken balance of phosphate repulsion 
at the opposite DNA sides that forces it to bend^l. If the 
A-tract curvature had an electrostatic origin, the phos- 
phate "hole" at the nick could perturb it by adding a 
local bend in a direction that should rotate as the "hole" 
is moved along the DNA chain. Like that the integral 
curvature would be deviated and partially compensated 
in a way compatible with the results shown in Fig. |2 

To check this possibility, similar experiments were car- 
ried out for a series of analogous phosphorylated nicks. 
Their relative PAGE mobilities are compared with other 
results in Fig. The mobility coefficients Q m used in 
this plate were computed as follows. The Q m value of 
the S fragment (see Table HJ is always the largest in the 
gel and it is arbitrarily assumed to equal 100. The Q m 
of fragment B is assumed to equal 0. For any given band 
its Q TO is promotional to the distance from the fastest S 
band in the same gel and it is estimated relative to the 
distance between B and S. If a nicked fragment migrates 
slower than the B-band in the same gel its Q m is neg- 
ative. In contrast, if it is faster, the corresponding Q m 
is positive below 100. Fig. displays thus obtained Q m 
values for all DNA fragments used in the present study. 

With free 5'-phosphates added to the nick sites, the 
phosphate "holes" are quenched or even inverted, which 
should drastically change the observed curvature modu- 
lations. The results in Fig. 03 however, show the oppo- 
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FIG. 3: Quantitative comparison of relative mobilities of 
DNA fragments listed in Table The method of calcula- 
tion of the mobility coefficients Q m is explained in the text. 
The reference mobilities of fragments B and S are indicated 
by horizontal dashed lines. The error bars indicate the esti- 
mate of experimental errors obtained from variation of Q m in 
repeated runs. 

site. The original skewed-sinusoidal dependence is persis- 
tent. At the same time, the phosphorylated nicks exhibit 
uniformly reduced mobilities, which agrees with earlier 
report ofiftifi 1 and should be attributed to the increased 
strand fraying in nicks. The additional phosphate charge 
increases the inter-strand electrostatic repulsion at sin- 
gle stranded breaks, which should enhance the fraying 
effect as discussed in detail elsewhere^. Whatever the 
reason of the uniformly reduced mobility, however, Fig. 
El clearly shows that the above hypothesis of electrostatic 
curvature compensation can be ruled out. 

This figure also shows the profile of mobilities of the 
nBo series in the presence of Mg 2+ ions. The Mg 2+ ions 
are long known to affect the DNA curvature by increasing 
or reducing it depending upon the sequenc e) 84 ' 85 . In our 
case the curvature was increased as judged from the ab- 
solute PAGE separation of fragments B and S. It is seen 
in Fig. that the profile of the nick mobilities in the 
presence of Mg 2+ ions remains qualitatively similar, but 
some nicks migrated slower than the mother B fragment. 
Interpretation of this effect depends upon the mechanism 
by which Mg 2+ increase the A-tract curvature. If the 
additional bend is caused by the electrostatic attraction 
due to Mg 2+ positioning between the phosphate groups 
in the narrowings of the major groove than the gain in 
the curvature may be larger when the DNA is nicked in 
the widening of the minor groove at the opposite side 
of the bend. This interpretation would agree with Fig. 
because the negative Q m are observed for nick posi- 
tions outside A-tracts in the probable widenings of the 
minor groove. It is also possible, however, that Mg 2+ 
ions somehow increase the nick fraying especially when 
the nick position approaches the poly (d A) .poly (dT) tails. 
In this case, the negative Q m values would have the same 
origin as for the phosphorylated nicks. 
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FIG. 4: PAGE analysis of comparative mobilities of double 
nicked DNA constructed from the random and A-tract repeat 
fragments. The contents of the lanes is as follows (see TableU 
and Methods for fragment codes). The lanes marked 'S,', 'B' 
and 'R' ("reference") contain the intact fragments S, B, and 
their mixture, respectively. Lanes marked by integer num- 
bers contain the corresponding single-nicked DNA fragments 
of the series marked below. Lanes marked 'D' contain DNA 
fragments carrying simultaneously the two nicks from the left 
and right neighbor lanes. 



The results in Fig. [21 can also be interpreted with the 
wedge model of DNA bending by assuming that nicks 
have a special structure with a wedge in a fixed local 
direction. The overall curvature is reduced when this 
wedge direction is opposite to that of the initial bend and 
increased in the opposite orientation. Similarly to the 
"phosphate hole" model above, skewed-sinusoidal modu- 
lations of curvature should occur as the nick site is trans- 
lated along DNA. Note that this wedge should have a 
very large angle since it compensates in nB-4o and nB+6o 
more than 30% of the curvature due to four phased A- 
tracts. A strong wedge like that should result in an in- 
creased curvature in nB+0o as well as in the nSo series in 
Fig. 01 which is not seen. This interpretation also would 
not agree with earlier studies of nick structures and gel 
mobilitie oftftifi 1 '^ . However, earlier PAGE studies were 
carried out at higher temperatures where the isotropic 
flexibility of nicks could, in principle, dominate all other 
effects. Also, the non-linear dependence of the gel mobil- 
ity upon the curvature may be such that the acceleration 
effect of the nick wedge in Fig. [2]is more pronounced than 
the expected retardation for nB+0o and the nSo series. 
Therefore, in order to check the above explanations, dou- 
ble nicked A-tract fragments were constructed, with the 
two sites separated by an odd number of helical half-turns 
so that the hypothetical wedges had opposite directions 
and should have compensated one another. The PAGE 
mobilities of double nicked fragments are compared in 
Fig. Q] These data evidently disagree with the forego- 
ing wedge interpretation of mobility modulations. For 
both nB-4o and nB+6o, the addition of a second nick at 
1.5 helical turn neither cancel nor even reduce the corre- 
sponding band shifts with respect to Fig. [2 As expected, 
analogous constructs for fragment S all show PAGE mo- 
bilities similar to that of the intact duplex. 

The experimental results presented above evidence 
that single stranded breaks really reduce the overall cur- 
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TABLE II: Structural parameters of standard and computed 
DNA conformations. The top row features the canonical B- 
DNA values corresponding to the fiber crystal structure—. 
The RMSD values were computed with respect to this con- 
formation. The results are shown for the last 1 ns averaged 
structures from the corresponding trajectories. The data for 
B and S trajectories are taken from our previous report^-. 
The helicoidals are the sequence averaged values computed 
with program Curves^. All distances are in angstroms and 
angles in degrees. 





Xdisp 


Inclin 


Rise 


Twist 


RMSD 


B-DNA 


-0.7 


-6.0 


3.4 


36.0 


0.0 


S 


-0.1 


-4.7 


3.5 


34.4 


4.1 


B 


-0.4 


-4.0 


3.5 


34.2 


6.8 


nB-4o 


-0.4 


-3.3 


3.5 


34.7 


6.0 


nB-2o 


-0.3 


-3.4 


3.5 


34.5 


3.2 


nB+Oo 


-0.1 


-5.0 


3.5 


34.5 


6.4 


nB+2o 


-0.7 


-4.4 


3.5 


34.1 


4.0 


nB+4o 


-0.6 


-4.5 


3.5 


34.1 


3.3 


nB+6o 


-0.6 


-4.3 


3.5 


34.3 


4.8 



vature of the A-tract repeat in remarkable qualitative 
agreement with CBT. Other theories of DNA bending 
cannot account for these data unless strong ad hoc as- 
sumptions are introduced like sequence-and-position de- 
pendent nick wedges, for example. 



TABLE III: Local inter-base parameters at strand breaks. 
The nicked steps only were analyzed with Curves^ in the last 
five Ins averaged structures. The same steps were analyzed 
in their intact states by similarly using the B trajectory of the 
mother A-tract fragment. Each fragment is characterized by 
two lines showing the averages and standard deviations for 
the given base step in its nicked (top) and intact (bottom) 
states. All distances are in angstroms and angles in degrees. 



Shift 


Slide 


Rise 


Tilt 


Roll 


Twist 








nB — 4o 






1.2±0.3 


-0.4±0.1 


3.3±0.1 


9.7±1.7 


-14. 6± 4.9 


40.5±1.3 


0.9±0.8 


-0.7±0.2 


3.2±0.2 


7.8±2.7 


-4.1± 5.1 


37.5±5.8 








nB - 2o 






0.5±0.2 


2.6±0.1 


3.2±0.1 


12.7±1.4 


4.2± 3.0 


7.8±2.9 


1.2±0.9 


0.1±0.2 


3.2±0.2 


10.3±3.4 


5.8± 2.7 


37.5±2.2 








nB + 0o 






-1.0+0.1 


1.2±0.2 


3.6±0.2 


-11.1±6.5 


1.8± 2.4 


52.2±2.8 


0.4±0.2 


-1.2±0.2 


3.5±0.1 


0.7±1.0 


-1.8± 2.7 


32.6±1.0 








nB + 2o 






-0.9±0.1 


-0.1±0.6 


2.8±0.1 


14.0±2.7 


-2.3± 7.9 


20.7±2.7 


0.1±0.2 


-1.8±0.1 


3.8±0.1 


-2.6±2.6 


-3.8± 1.5 


32. Oil. 6 








nB + 4o 






-0.6±0.2 


-1.3±0.2 


3.3±0.2 


5.3±2.8 


-10.8±10.1 


25.6±1.7 


-0.2±0.1 


-1.5±0.2 


3.4±0.1 


3.6±2.1 


-10. 3± 2.8 


30.4±0.8 








nB + 6o 






0.8±0.1 


-0.4±0.3 


3.4±0.2 


9.6±2.0 


-23. 9± 5.1 


39.0±1.7 


0.5±0.2 


-1.1±0.2 


3.4±0.2 


3.7±3.0 


-8.3± 1.2 


37.0±2.1 



C. DNA Curvature in Simulations 

As regards the overall stability of dynamics and its 
closeness to the B form, all MD trajectories computed in 
the course of this studies were similar to our earlier sim- 
ulations under similar conditions£2iitiii£i2tii2i. Table |n] 
shows parameters of the final Ins-average conformations. 
They all have remarkably similar average helicoidals cor- 
responding to a typical B-DNA. For example, the average 
helical twist estimated from the best-fit B-DNA experi- 
mental values* gives 34.0 ± 0.2° and 33.8 ± 0.2° for the 
A-tract fragment and the randomized sequence, respec- 
tively. The nicks did not cause readily visible perturba- 
tions and it was commonly difficult to distinguish them 
in snapshots. Statistical analysis of the local inter-base 
parameters at nick positions is summarized in Table ITTT1 
As expected, the broken DNA strand exhibits somewhat 
increased flexibility at the nick site. The most noticeable 
are much stronger deviations of Twist form its canonical 
values. A similar feature is smaller, but distinguishable 
for Tilt and Roll. On the other hand, the nanosecond 
time scale fluctuations characterized by the standard de- 
viations are not very different with and without the nick. 
The last observation suggests that, as long as the bases 
at the nicked step remain in the stack, their rapid mo- 
tions are restrained to the same degree as in the intact 
structure, and that the above strong deviations of some 
averages are due to slower global motions that occur in 
much longer time scales. Apart from these strong fluctu- 



ations, however, we did not observe any repetitive static 
perturbations that might be attributed to strand break- 
ing. Moreover, strongly non-canonical helical parameters 
are sometimes encountered in intact steps as well. As it 
was for dynamics of the mother A-tract DNA fragment, 
the rms deviations from the canonical structures seen in 
Table |n] are largely due to variation of the overall curva- 
ture. 

Fig. compares the overall DNA curvature in the 
last five Ins averaged conformations from all trajectories 
mentioned in Table ITfl It is seen that the mother A-tract 
fragment (B) reached the largest bending amplitudes non 
accessible to other fragments. The most stable and signif- 
icant curvature among the nicks is obtained for nB-4o and 
nB+0o. Detailed in Fig. is the dynamics of the bend 
direction. Its short time scale fluctuations are smaller 
when bending is strong, therefore, Fig. El also gives an 
estimate of the bending amplitude. Again we note that 
the mother A-tract fragment (B) is distinguished by the 
rapid establishment of curvature. However, convergence 
to a statically bent state is clear for nB-4o and nB+0o 
as well. Following to our previous approach2SiiSii& an 
attempt has been made to confirm the attracting prop- 
erty of the bent state for nB+0o by running dynamics 
starting from the canonical A-DNA conformation. The 
corresponding trajectory was continued to 17 ns and it 
showed convergence of the same bending direction as in 
Fig. El but with a somewhat smaller amplitude. 

Fig. Eland Fig. EOshow that our MD simulations failed 
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FIG. 5: (a) Geometric constructions used for evaluating the 
DNA bending. The amplitude of bending is measured by the 
angle between the ends of the curved helical axis produced by 
the Curves algorithm—. The two coordinate frames shown 
are the global Cartesian coordinates (OXYZ), and the local 
frame constructed in the middle point according to the Cam- 
bridge convention (O'JKL)— . The curve is rotated with two 
its ends fixed at the Z-axis to put the middle point in plane 
XOZ. The bending direction is measured by angle (fi between 
this plane and vector J of the local frame. By definition, this 
vector points to the major DNA groove along the short axis 
of the reference base pais^. Consequently, the zero <p value 
corresponds to the overall bend towards the minor groove in 
the middle of the DNA fragment as it was in the very first 
analyses of local DNA curvature^. 

(b) Polar plot of bending angle (A) versus direction (phi). 
The results are shown for the last five 1 ns averaged struc- 
tures from each trajectory, with the curvature measured as 
explained in plate (a). 



to reproduce a regular variation of the bending magni- 
tude with the nick position as seen in the PAGE exper- 
iments above. Fragments nB-4o and nB+0o exhibit no- 
ticeable curvature comparable, but smaller than that of 
the intact duplex. According to experiment, the curva- 
ture in nB+0o is similar to that in the mother fragment, 
but for nB-4o it should be definitely smaller. Fragment 
nB+6o exhibited less stable bending in a similar direction 
whereas for nB-2o, nB+2o, and nB+4o the 20 ns duration 
of trajectories apparently was not sufficient for conver- 
gence. At the same time, our calculations confirm that, 
in qualitative agreement with experiment, single strand 
breaking interferes with bending and generally reduces 
it. It is remarkable in Fig. EJthat strong bends beyond 
40° were observed only in directions roughly correspond- 
ing to that in the mother fragment, which confirms the 
strong attracting property of this bent state. Taking into 
account the unclear physical nature of the phenomenon 
being studied and the complexity of the model system we 
consider the partial qualitative agreement achieved here 
as very satisfactory. 

In spite of this only limited agreement with experi- 
ment, we tried to use the ensemble of computed curved 
DNA conformations to clarify the origin of the effect pro- 
duced by nicks upon the A-tract curvature. To this end, 
for each of the fragments B, nB-4o, nB+0o, and nB+6o a 
nanosecond of dynamics have been selected for a more de- 
tailed analysis. The corresponding Ins average structures 
are marked by arrows in Fig. They all correspond to 
the largest degree of curvature reached during these tra- 
jectories. For B and nB+6o these are the last nanosec- 
ond structures. In contrast, for nB-4o and nB+0o the 
17th and the 19th nanoseconds were used, respectively. 
According to Fig. Q[J these bends are more or less conver- 
gent suggesting that the corresponding dynamics should 
be similar and subtle differences caused by nicks might 
be detectable. As seen in Fig. for fragments B, nB-4o, 
and nB+6o the averaged conformations are rather close 
to one another. The curved DNA axis is almost planar 
with two or three zones of increased bending found be- 
tween the A-tracts. In contrast, for nB+0o the curved 
axis is strongly non-planar due to significant deviation 
of bend direction in its central zone from the other two. 
As a result, for B, nB-4o, and nB+6o the overall bend is 
close to the sum of local bends, but for nB+Oo it is much 
smaller. The minor groove profiles averaged over the cor- 
responding 1 ns periods are compared in Fig. 03 Again 
one can note a similarity between B, nB-4o, and nB+6o 
versus a more complex profile for nB+Oo. The widenings 
of the minor groove are found between A-tracts and they 
roughly correspond to the zones of increased bending in 
Fig. [3 According to Fig. [3 and Fig. |H1 a reduced overall 
bending amplitude in nB-4o and nB+6o with respect to 
the mother A-tract fragment is attributable to the cen- 
tral zone where the widening is also less significant. In 
contrast, the smaller curvature in nB+0o is mainly due 
to the lost phasing between the local bends. 

Relaxation of A-tract curvature by single stranded 
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breaks has been originally proposed as a conceptual test 
of CBT22.. The results reported here qualitatively agree 
with its predictions and strongly support this model ver- 
sus its alternatives. At the same time, the nick position 
dependence as revealed by our PAGE tests is inverted 
with respect to to the original guessS. We thought that 
the backbone compression is maximal in the widenings 
of the minor groove, therefore, nicks between A-tracts 
should relax bending stronger than nicks inside them. 
In contrast, the experiments reported above suggest that 
the backbone compression should increase inside A-tracts 
and reduce outside them. Because A-tracts are found 
at the internal edge of the curved double helix the last 
suggestion agrees with the simple physical intuition that 
says that the surface of a curved cylinder is compressed 
at its inner edge and stretched at the opposite side. In 
order to understand which backbone component can be 
responsible for the compression/frustration relationship 
postulated by CBT we tried to measure the specific back- 
bone length by using pairs of similar atoms in neighboring 
residues. 

For all atom pairs the resulting profiles of the backbone 
length appear very noisy and regular modulations can be 
seen only after averaging with a sliding window. Inspec- 
tion of all atom pairs allows one to divide them in two 
groups. For atoms of phosphate groups and their neigh- 
bors, the measured backbone length usually is larger at 
the inner side of the bend and smaller at the opposite 
side, which may be viewed as a trend towards the A- 
from in the widenings of the major groove. However, 
these modulations are very weak, in agreement with the 
conventional view of the B-DNA backbone as being non- 
compressibl o 90 i 91 . In contrast, for a group of backbone 
atoms including 04' and its close neighbors, strong regu- 
lar oscillations of the backbone length could be revealed, 
with their phases being opposite to those in the first 
group. The average backbone profiles produced by inter- 
04' distances are shown in Fig. They were computed 
for the 1 ns intervals of trajectories marked in Fig. It 
is seen that the backbone length correlates with the mi- 
nor groove profiles in Fig. |H1 but does not repeat them. 
Thus measured specific backbone length reaches its lo- 
cal minima at the inner edge of the bent DNA cylinder 
and corresponds exactly to the nick positions in nB-4o 
and nB+6o. Similar oscillations are observed in X-ray 
structures of bent DNA available in NDB22. (AKM, un- 
published). It may seem surprising that profiles in Fig. 
El exhibit no perturbations at strand breaks, but this re- 
sults from smoothing. In the non-smoothed plots the 
single stranded breaks in nB+6o and nB+Oo were dis- 
tinguished by very small and very large distances, re- 
spectively, but they have been compensated by strong 
opposite deviations in the neighboring steps. 

The foregoing analysis demonstrates that the DNA 
backbone behaves as a complex mechanical system rather 
than a chain of strings or rigid rods, and that different 
its components compress and stretch in different zones of 
the curved DNA. Nevertheless, the correlations revealed 



in Fig. Elsuggest that the backbone compression, if it re- 
ally exists, should mainly affect sugar-sugar interactions 
and that they can be at the origin of the intrinsic back- 
bone frustration. Even though in B-DNA consecutive 
04' atoms do not interact, the neighboring sugar rings 
make a number of direct contacts non-mediated by the 
solvent. 



IV. CONCLUDING DISCUSSION 

The results presented here provide new insights in 
the putative physical mechanism of intrinsic DNA cur- 
vature. This intriguing phenomenon seems to be exhaus- 
tively studied, but it still attracts great interest. The 
A-tract repeats and single stranded breaks both have 
been earlier recognized as elements specifically involved 
in DNA bending. To our best knowledge, however, no- 
body checked what could be their cumulative effect. Here 
it is found that nicks relax the intrinsic curvature induced 
by A-tract repeats in conditions where they do not affect 
the random sequence DNA. The relaxation effect depends 
upon the nick position with respect to the bend. It is very 
significant when the nick is at the inner edge of the bent 
DNA inside A-tracts, and is gradually reduced essentially 
to zero when the nick is moved to the outer surface of the 
bend. This behavior is in qualitative agreement with the 
CBT and it suggests that the backbone compression in 
the curved DNA is larger in the minor groove narrowings 
at the inner surface. 

The MD simulations confirm that single stranded 
breaks interfere with, and tend to relax the A-tract in- 
duced curvature. For about a half of the DNA frag- 
ments studied, the duration of trajectories provided re- 
producible convergence to statically bent states in good 
agreement with experiment as well as earlier simulations. 
Inspection of the computed conformations reveals strong 
regular modulations of the local backbone length as mea- 
sured by distances between certain sugar atoms, with 
compression reaching its maximum at the inner edge of 
the bend where nicks produce the strongest relaxation in 
experiments. These results provide additional support to 
the CBT and suggest that the frustration in the B-DNA 
backbone may result from interactions between consecu- 
tive sugar rings. 
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FIG. 6: The time evolution of the bending direction as mea- 
sured by the <p angle in Fig. (in degrees). Here the axis 
of the curved double helix was computed as the best fit com- 
mon axis of coaxial cylindrical surfaces passing through sugar 
atoms, which gives solutions close to those produced by the 
Curves algorithm^. The traces have been smoothed by aver- 
aging with a window of 60 ps and are displayed with a time 
step of 40 ps. 



FIG. 7: Schematic drawings of the of 1 ns averaged struc- 
tures marked by arrows in Fig. The pictures were pro- 
duced with Curves^. The backbone breaks are distinguished 
by three consecutive gaps in the ribbons, but only the central 
phosphate is in fact absent. Two orthogonal views are shown. 
The upper and lower views correspond to the negative direc- 
tion of the OY axis and the positive direction of the OX axis, 
respectively, as shown in in Fig. Ht. 
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FIG. 9: Backbone length profiles averaged over the 1 ns tra- 
jectory intervals marked in Fig. The plots were generated 
by using inter-04' distances. The central traces represent the 
average with rms fluctuations shown as error bars. The upper 
and lower solid traces show the maximal and minimal values, 
respectively. All profiles were smoothed with a sliding window 
of 3 base pair steps. The length is given in angstroms, with 
the corresponding canonical A- and B-DNA levels indicated 
by the horizontal dashed lines. The vertical grey bars mark 
nick positions. 



FIG. 8: The profiles of the minor groove profile averaged 
over the 1 ns trajectory intervals marked in Fig. [£] The 
central traces represent the average groove width with rms 
fluctuations shown as error bars. The upper and lower solid 
traces show the maximal and minimal values, respectively. 
The groove width is evaluated by using space traces of C5' 
atoms—. Its value is given in angstroms, with the correspond- 
ing canonical A- and B-DNA levels indicated by the horizontal 
dashed lines. The vertical grey bars mark nick positions. 
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